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Abstract 

In this note, we study the phase transitions arising in a modified Smolu- 
chowski equation on the sphere with dipolar potential. This equation models 
the competition between alignment and diffusion, and the modification con- 
sists in taking the strength of alignment and the intensity of the diffusion as 
functions of the order parameter. 

We characterize the stable and unstable equilibrium states. For stable 
equilibria, we provide the exponential rate of convergence. We detail special 
cases, giving rise to second order and first order phase transitions, respectively. 
We study the hysteresis diagram, and provide numerical illustrations of this 
phenomena. 

1 Introduction 

In this short note, we study the following modified Smoluchowski equation (also 
called Fokker-Planck equation), for an orientation distribution f{u},t) defined for a 



1 



time t ^ and a direction u; G S (the unit sphere S of M") as follows: 

dtf = -iy{\Jf\)V^ ■ ifVUoo ■ m + T{\Jj\)A^f =: Qif), (1) 

fi/ = ni' •^/(^)= / (2) 

where A^^, V^j-, and V^; are the Laplace-Beltrami, divergence, and gradient oper- 
ators on the sphere. The vector Jf G is the first moment associated to / (the 
measure on the sphere is the uniform measure such that j^dv = 1), and ^2/ G S 
represents the mean direction of the distribution /. 

The first term of the right-hand side of ([T]) is an alignment term towards Qf, 
and the function u represents the strength of this alignment. The function r is the 
intensity of the diffusion on the sphere. 

When r is a constant and i^(|J/|) = we can recover the standard Smolu- 
chowski equation on the sphere, with dipolar potential [TT]. Indeed, the dipolar 
potential is given by K{u,u) = —u ■ u, and the equation can be recast as: 

dtf = V^- (^fV^{J^K{uj-oo)f{oo,t)du)^+rA^f. (3) 

Other classical kernels [HlIlllIIS] for the study of semi-dilute and concentrated 
suspensions of polymers are the so-called Maier-Saupe potential K{uj, u) = — (w-o;)^, 
or the original Onsager potential K{u, u) = \u}xu\. In these cases, there are a lot of 
studies regarding the phase transition phenomenon for equilibrium states Ptl^l^fTUj 
fT2 | [T3 |[T6] - [20] ■ and in particular a hysteresis phenomenon occurs for the Maier-Saupe 
potential in dimension 3. The case of the dipolar potential has also been studied 
precisely [11] , with an analysis of the rates of convergence to the equilibrium as time 
goes to infinity. In this case, there exists a so-called continuous phase transition 
for a critical threshold t^. when t < Tc the solution converges exponentially fast 
to a non-isotropic equilibrium; when r > Tc it converges exponentially fast to the 
uniform distribution. At the critical value r = Tc, the solution converges to the 
uniform distribution at a rate 

Here we study the modifications arising when i/ and r depend on \Jf\, which 
can be motivated by some biological modeling (see [7j and references therein). In 
that case, we cannot use r as a bifurcation parameter anymore. Instead, we will use 
the initial mass p (a conserved quantity) as the key parameter to study the phase 
transition. We will assume that: 

Hypothesis 1.1. 

(i) The functions v and r are , with z/(0) = and r > 0. 

(ii) The function \ J\ i— t- h{\J\) = is an increasing function. We denote by a 

its inverse, i.e. 

K = h{\J\) ^ \ J\ = a{K). 

The first part of this hypothesis implies that we do not have any singularity 
of Q in (II]) as I J/-| — 7- 0: if | J/| = 0, we simply have Q{f) = r(0)Ai^/. The second 
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part is made for the sake of simplicity. It leaves enough flexibility to to reveal key 
behaviors in terms of phase transitions. It would be easy to remove it at the price of 
an increased technicality. Additionally, it means that when / is more concentrated 
in the direction of Qf, the relative strength of the alignment force compared to 
diffusion is increased as well. This can be biologically motivated by the existence of 
some social reinforcement mechanism. 

We will see that we can observe a wealth of phenomena, including hysteresis. 
The purpose of this note is to summarize the analytical results, as well as some 
numerical simulations which illustrate this phenomena. All the proofs are detailed 
in mj, where ([T]) arises as the spatially homogeneous version of a space-dependent 
kinetic equation, obtained as the mean-field limit of a self-propelled particle sys- 
tem interacting through alignment. This spatially homogeneous study is crucial to 
determine the macroscopic behavior of this space-dependent kinetic equation. 

2 General study 

2.1 Existence and uniqueness 

We first state results about existence, uniqueness, positivity and regularity of the 
solutions of (II]). Under hypothesis 11.11 we have the following 

Theorem 1. Given an initial finite non-negative measure /o in iJ*(S), there exists 
a unique weak solution / o/ ([T]) such that /(O) = /q. This solution is global in time. 

Moreover, f G C°°((0, +oo) x E>), with f{t,u) > for all positive t, and we have 
the following instantaneous regularity and uniform boundedness estimates (for m G 
N, the constant C being independent of f), for all t > 0: 



For later usage, we define as an anti-derivative of h: = h{\J\). In 

this case, the dynamics of corresponds to the gradient flow of the following free 
energy functional: 





Indeed, if we define the dissipation term X'(/) by 



V{f)=T{\Jf\) ! /|V.(ln/-M|J;|)a;-fi;)|Ma;, 



we get the following conservation relation: 



d_ 



■Hf) 



v{f) <: 0. 



(4) 
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2.2 Equilibria 

We now define the von Mises distribution which provides the general shape of the 
non-isotropic equilibria of Q. 

Definition 2.1. The von Mises distribution of orientation G S and concentration 
parameter k, is given by: 

The order parameter c{k,) is defined by the relation 
and has expression: 

rcos^e'^^°^^sin"-2 0d^ 

c(k) = ^ . 6 

^ ' /o"e«<=°^^sin"~2^d^ ^ ' 

The concentration parameter c(k) defines a one-to-one correspondence k G [0, oo) H- 
c(k) G [0, 1). The case k = c(k) = corresponds to the uniform distribution, while 
when K is large (or c{k) is close to 1), the von Mises distribution is closed to a Dirac 
delta mass at the point Q. 

Some comments are necessary about the interval of definition of a. First note 
that, under hypothesis 11.1^ h is defined from [0,+oo), with values in an inter- 
val [0,Kmax), where we may have Kmax = +oo. So a is an increasing function 
from [0, Kmax) onto M+. Moreover, for later usage, we can define 

ro = r(0)>0, and Pe=lim^ = ^ (7) 

where Pc > may be equal to +oo, and where we recall that n denotes the dimension. 
The equilibria are given by the following proposition: 

Proposition 2.1. The following statements are equivalent: 

(i) f G C2(§) and Qif) = 0. 

(ii) f G C\S) andVW) = 0. 

(Hi) There exists p ^ and f2 G S such that f = pM^Q, where k ^ satisfies the 
compatibility equation: 

(t(k) = pc{k). (8) 

Let us first remark that the uniform distribution, corresponding to k = is 
always an equilibrium. Indeed, we have c(0) = a{fS) = and is satisfied. How- 
ever, Proposition 12.11 does not provide any information about the number of the 
non-isotropic equilibria. Indeed, equation can be recast into: 

4^ = i. (9) 

a{K) p 
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which is vahd as long as a 7^ 0. We know that cr is an increasing unbounded function 
from its interval of definition [O^Hmax) onto [0,+oo), and thanks to hypothesis 11.11 
and to (171), we know that ~ as k — (if pc < +00). So since c(k) ~ 
as K — > (see for instance [H]), we have the two following results (also valid in the 
case Pc = +00): 



c(k) 1 c(k) 

— -— — )■ — as K — and — -— — as k — )■ 



a{K) Pc a[K) 
We deduce that this function reaches its maximum, and we define 



(10) 



p* 



mm 



c{k) 



(11) 



For p < p^, the only solution to the compatibility condition is k = 0, and the 
only equilibrium is the uniform distribution f = p. Except from these facts, we 
have no further direct information of this function k i— )■ c{K)/a{K), since c and a 
are both increasing. Figure [T] depicts some examples of the possible shapes of the 
function n H- c(fi:)/cr(fi:). 
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Figure 1: The green, blue, red and purple curves correspond to various possible 
profiles for the function k i— )■ 



We see that depending on the value of p, the number of families of non-isotropic 
equilibria, given by the number of positive solutions of the equation ([9]), can be zero, 
one, two or even more. We now turn to the study of the stability of these equilibria, 
through the study of the rates of convergence. 
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2.3 Rates of convergence to equilibrium 

The main tool to prove convergence of the solution to a steady state is LaSalle's 
principle, that we recall here (the proof follows exactly the lines of [H]). By the 
conservation relation (jl]), we know that the free energy T is decreasing in time 
(and bounded from below since | J| is bounded). LaSalle's principle states that the 
limiting value of T corresponds to an w-limit set of equilibria: 

Proposition 2.2. LaSalle's invariance principle. 

Let /o be a positive measure on the sphere S. We denote by Too the limit 
of J^{f{t)) as t ^ oo, where f is the solution to the modified Smoluchowski equa- 
tion ([T]) with initial condition /q. 

Then the set £^ = {f e C°°{S) s.t. V{f) = and F{f) = T^} is not empty. 

Furthermore f{t) converges in any norm to this set of equilibria (in the 
following sense): 

lim inf \\f{t)-g\\Hs = Q. 

Since we know the types of equilibria, we can refine this principle to adapt it to 
our problem: 

Proposition 2.3. Let /q be a positive measure on the sphere §, with initial mass p. 

If no open interval is included in the set {n, pc{k) = a{K)}, then there exists a 
solution Koo to the compatibility solution such that we have: 



lim \Jf{t)\ = pc(Koo) 
\fs e R, lim ||/(t) - pM,^nfit)\\H^ = 0. 



This last proposition helps us to characterize the w-limit set by studying the 
single compatibility equation ([H]). 

When K = is the unique solution, then this gives us that / converges to the 
uniform distribution. Otherwise, two cases are possible, either = 0, and / 
converges to the uniform distribution, or 7^ 0, and the only unknown behavior 
is the one of ^f{t)- If we are able to prove that it converges to fioo £ S, then / 
converges to a fixed non-isotropic steady-state pM^^Q^. 

However, Proposition 12.31 does not give information about quantitative rates of 
convergence of | J/| to pc{Koo), and of ||/(t) — pM,^^fi^.(^t)\\H'' to 0, as t — 00. So we 
now turn to the study of the behavior of the difference between the solution / and 
a target equilibrium pMf,^a^(^t)- 

This study consists in two types of expansion. If we expand the solution around 
the uniform equilibrium, some simple energy estimates give us exponential con- 
vergence when p < pc- But when we expand the solution around a non-isotropic 
equilibrium pM^^Q^(^t), we see that the condition of stability is related to the mono- 
tonicity of the function k i— c{K,)/a{K). Hence we can see directly on the graph 
of this function (see examples on Figure [1]) both the number of family of equilibria 
and their stability: if the function is decreasing, the family is stable. By con- 
trast it is unstable when the function is increasing. When the difference between / 
and pMK^nj(t) converges exponentially fast to (on the stable branch), we are able 
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to control the displacement of ^f{t), which gives convergence to Qoo ^ S. We then 
have convergence of / to a given equilibriumpMK^Qoo- 

All these results are summarized in the following two theorems. In what follows, 
we say that a constant is a universal constant when it does not depend on the 
initial condition /o (that is to say, it depends only on p, n and the coefficients of 
the equation u and r, and on the exponent s of the Sobolev space in which the 
result is stated). 

Theorem 2. We have the following instability and exponential stability results 
around the uniform equilibrium: 

• Suppose that p < pc- We define 

A = (n-l)ro(l-^)>0. (12) 

pc 

There exists a universal constant C, such that if \\fo — pWh'^ < ^, then for 
all t 0, we have 

l-fll/o-p|k« 



\\f{t)-p\\H^ ^ e 



• If p > Pc, and if Jf^ ^ 0, then we cannot have k^q = in Proposition }^.!^ - the 
solution cannot converge to the uniform equilibrium. 

To study the stability around a non-isotropic equilibrium, we fix p, and we de- 
note by K a positive solution to the compatibility equation (we will not write the 
dependence of c and a on k when there is no possible confusion). We denote by J-'^ 
the value of J-'(pMkq) (independent of G §). 

Theorem 3. We have the following instability and exponential stability results when 
starting close to a non-isotropic equilibrium: 

• Suppose (f )'(k) > 0. For all s > there exist universal constants 6 > 
and C > 0, such that for any initial condition /o satisfying \\fo — pMi^^W^^ < S 
for some f2 G S, there exists Q^o G S such that 



-xt 



\\f-pM^nJ\H^ ^CWfo-pM^njH^e- 
where the rate is given by 

A = ^A.A'. (13) 

a' c 

The constant is the best constant for the following weighted Poincare in- 
equality (see the appendix of for more details on this constant, which does 
not depend on fl): 

I I V^(7|' M^n du^Aj I g^M^n dco - { [ g M^n dw)'! • (14) 
JS '-Js Js 

Suppose (^)'(/t) < 0. Then any equilibrium of the form pM^Q is unstable, 
in the following sense: in any neighborhood of pM^^Q, there exists an initial 
condition /q such that J^{fo) < Tn- Consequently, in that case, we cannot 
have Proposition \2.!A 
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3 Second order phase transition 



Let us now focus on the case where we always have (^)' > for all k > (see 
for example the lowest two curves of Figure [T]). In this case, the compatibility 
equation Q has a unique positive solution for p > p^. With the results of the 
previous subsection about stability and rates of convergence, we obtain the behavior 
of the solution for any initial condition /o with initial mass p. 

• If p < Pc then the solution converges exponentially fast towards the uniform 
distribution f^o = p- 

• If p = Pc the solution converges to the uniform distribution. 

• If p > Pc and Jfy 7^ 0, then there floo such that / converges exponentially fast 
to the von Mises distribution f^o = pMi^^^ , where k > is the unique positive 
solution to the equation pc(k) = cr(K). 

The special case where J/q = leads to the heat equation dtf = r^^^f . Its solution 
converges exponentially fast to the uniform distribution, but this solution is not 
stable under small perturbation of the initial condition. Let us remark that for some 
particular choice of the coefficients, as in [H], it is also possible to get an algebraic 
rate of convergence in the second case p = pc- For example when (j{k) = n, we 
have 11/ — pII ^ Ct~2 for t sufficiently large. 

So we can describe the phase transition phenomena by studying the order param- 
eter of the asymptotic equilibrium c = as a function of the initial density p. 

We have c(p) = if p ^ pc, and c is a positive continuous increasing function 
for p > Pc. In the common situation where ~ = ~ — '^^'^ + o{k7) when k — )■ 0, it is 
easy to see, since c{k) ~ when k — )■ 0, that we have 

c(p) ~ a(p - pc)^, as p 4 Pc. (15) 

Since ^ is Lipschitz, we always have /3 ^ 1. So the first derivative of c is discontinuous 
at p = Pc- This is the case of a second order phase transition (also called continuous 
phase transition). The critical exponent /3 can take arbitrary values in (0, 1], as can 
be seen by taking h{\J\) such that a^n) = c(k)(1 + k?). 

In general, we have the following practical criterion, which ensures a second order 
phase transition. 

Lemma 1. // is a non-increasing function of \J\, then we have (^)' > for 
all n > 0. In this case, the critical exponent (5 in (fT^ . if it exists, can only take 
values in [|, 1] ■ 

4 Hysteresis 

4.1 Typical example 

We now turn to a specific example, where all the features presented in the stability 
study can be seen. We focus on the case where i^(| J|) = as in [11], but we now 
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take r(| J|) = 1/(1 + \ J\)- From the modeling point of view, this occurs in the Vicsek 
model with vectorial noise (also called extrinsic noise) [HE]- 

In this case, we have /i(|J|) = |J| + |Jp, so the assumptions of Lemma [Hare not 
fulfilled, and the function a is given by a^n) = |(-\/l + — 1). 

Expanding ^ when k, is large or k, is close to 0, we get 

c _ + ^K + 0(/«2) asft^O, 
a |--^ + 0(k~^) as ft — )■ cxd. 

Consequently, there exist more than one family of non-isotropic equilibria when p is 
close to pc = n (and p > pc). 

The function k i— )■ can be computed numerically. The results are displayed 
in Figure |2] in dimensions n = 2 and n = 3. 

1 I 1 1 1 1 





0.2 - 



Figure 2: The function k H- -4^, in dimensions 2 and 3. 



We observe the following features: 

• There exists a unique critical point ft* for the function -, corresponding to 
its global maximum y- (in dimension 2, we obtain numerically p* ~ 1.3726 
and ft* ~ 1.2619, in dimension 3 we get p* ~ 1.8602 and ft* ~ 1.9014). 

• The function ^ is strictly increasing in [0, ft*) and strictly decreasing on (ft*, oo). 

From these properties, it follows that the solution associated to an initial condi- 
tion /o with mass p can exhibit different types of behavior, depending on the three 
following regimes for p. 
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• If p < p*, the solution converges exponentially fast to the uniform equilib- 
rium /oo = p. 

• If p* < p < n, there are two families of stable solutions: either the uniform 
equilibrium / = p or the von Mises distributions of the form pM^j^, for f2 G S 
where k is the unique solution with /t > of the compatibility equation (jH]). 
If /o is sufficiently close to one of these equilibria, there is exponential conver- 
gence to an equilibrium of the same family. 

The von Mises distributions of the other family (corresponding to solution 
of (IH]) such that < k < k*) are unstable in the sense given in Theorem [3l 

• If p > 77, and J/p 7^ 0, then there exists Vtao G S such that / converges exponen- 
tially fast to the von Mises distribution pMi^^^ , where n is the unique positive 
solution to the compatibility equation pc(k) = a{K,). 

At the critical point p = p^, the uniform equilibrium is stable (and for any initial 
condition sufficiently close to it, the solution converges exponentially fast to it), but 
the stability of the family of von Mises distribution {p*M^^^, G S} is unknown.. 

At the critical point p = n, the family of von Mises distribution {riM^^Q, f2 G S} 
is stable, where Kc is the unique positive solution of (JHl). For any initial condition 
sufficiently close to uMi^^q for some G S, there exists Qoo such that the solution 
converges exponentially fast to tt-M^^q^. However, in this case, the stability of the 
uniform distribution / = ra is unknown. 

As previously, in the special case Jf^ = 0, the equation reduces to the heat 
equation and the solution converges to the uniform equilibrium. 

Since c{k) is an increasing function of k, we can invert this relation k, h- c(k) 
into c h-> k(c) and express the density p = '^^'^^^^•'^ as a function of c. The result is 
depicted in Figure |3] for dimension 2 or 3. With this picture, we recover the phase 
diagram in a conventional way: the possible order parameters c for the different 
equilibria are given as functions of p. The dashed lines corresponds to branches of 
equilibria which are unstable. 

We can also obtain the corresponding diagrams for the free energy and the rates 
of convergences. For this particular example, the free energies J^{p) and J-'^ (we 
recall that they correspond respectively to the free energy of the uniform distribution 
and of a von Mises distribution pM^n for a positive solution k of the compatibility 
equation ([8]), including both stable and unstable branches) are given by 



J^(p) = plnp, 
J'^, = p In p + (p In M^q)m - -cr^ - -cr^ 

/I 2 
gKcosS^^ (k- a) + -an. 
6 3 

The plots of these functions are depicted in dimensions 2 and 3 are depicted on the 
left part of Figure HI Since the functions are very close in the figure for some range of 
interest, we depict the difference J-'^ — J-'(p) in a more appropriate scale, in the right 
part of Figure HI The dashed lines correspond to unstable branches of equilibria. 
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Density p 



Figure 3: Pliase diagram of the model with hysteresis, in dimensions 2 and 3. 




Density p p* pi p* Pi Density p 

Figure 4: Free energy levels of the different equilibria (left), and difference of free 
energies J-'^ — J^{p) (right), as functions of the density, in dimensions 2 and 3. 



We observe that the free energy of the unstable non-isotropic equilibria (in dashed 
line) is always above that of the uniform distribution. There exist pi G {p^,pc) 
and a corresponding solution k,i of the compatibility solution ([H]) (with ki > k*, 
corresponding to a stable family of non-isotropic equilibria) such that J-'^^ = J^{pjr). 
If p < Pi, the global minimizer of the free energy is the uniform distribution, while 
if p > pi, then the global minimum is reached for the family of stable von Mises 
equilibria. The physical relevance of this value is not clear though, as we will see in 
the numerical illustration of the next subsection. 

The rates of convergence to the stable equilibria, following Theorems |2] and [3l 
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are given by 



Ao = (n - 1)(1 - -), for p < Pc 
n 



1 1 n — 1 
K = — — A«(l - (- - c )(t(1 + 2a)), for p > p,, 

1 + 0" C K 

where Aq is the rate of convergence to the uniform distribution p, and A^ is the rate 
of convergence to the stable family of von Mises distributions pM^^, where k is the 
unique solution of the compatibility condition ([H]) such that k > n^,. Details for the 
numerical computation of the Poincare constant are given in the appendix of [S]. 
The computations in dimensions 2 and 3 are depicted in Figure [5l 



CD 



O 

CD 

cd 
Pi 



p* 2 3 
Density p 



Figure 5: Rates of convergence to both types of stable equilibria, as functions of the 
density p, in dimensions 2 and 3. 



4.2 Numerical illustrations of the hysteresis phenomenon 

In order to highlight the role of the density p as the key parameter for this phase 
transition, we introduce the probability density function f = j; and we get 

dtf = T{p\Jjr\)Aj- KpI^/I)V. ■ ifV^in^u)). (16) 

When p is constant, this equation is equivalent to ([T]). We now consider p as a 
parameter varying slowly with time (compared to the time scale of the convergence to 
equilibrium, see Figure E]). If this parameter starts from a value p < p*, and increases 
slowly, the only stable distribution is initially the uniform distribution / = 1, and it 
remains stable. So we expect that the solution stays close to it, until p reaches the 
critical value pc- For p > pc, the only stable equilibria are the von Mises distributions, 
and the solution converges to one of these equilibria. The order parameter defined 
as c(/) = I Jr|, then jumps from to Cc = c(kc). If then the density p is further 
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decreased slowly, the solution stays close to a von Mises distribution, and the order 
parameter slowly decreases, until p reaches p* back. For p < p^, the only stable 
equilibrium is the uniform distribution, and the concentration parameter jumps 
from c^: = c(k*) to 0. This is a hysteresis phenomenon: the concentration parameter 
describes an oriented loop called hysteresis loop. 

Let us now present some numerical simulations of the system (|TB|) in dimen- 
sion n = 2. We start with a initial condition which is a small perturbation of the 
uniform distribution, and we take p = 1.75 — 0.75cos(^t), with T = 500. We use a 
standard central finite different scheme (with 100 discretization points), implicit in 
time (with a time step of 0.01). The only problem with this approach is that the 
solution converges so strongly to the uniform distribution for p < p^, so after pass- 
ing pc, the linear rate of explosion for Jj is given by e^p^~^^^, and is very slow when p 
is close to Pc- So since Jj is initially very small when passing the threshold p = Pc 
we would have to wait extremely long in order to see the convergence to the stable 
von Mises distribution. To overcome this problem, we adding a threshold e and 
strengthen \Jjr\ when ||/ — l||oo^£, by 

f^f + max(0,£ - ||7- l||oo) 

We note that this transformation that we still have || / — 1 ||oo ^ £^ if it was the case 
before applying the transformation. 

Figure |6] depicts the result of a numerical simulation with a threshold e = 0.02. 
We clearly see this hysteresis cycle, which agrees very well with the theoretical 
diagram. The jumps at p = and p = Pc are closer to the theoretical jumps when T 
is very large. We were not able to see any numerical significance of the value pi (for 
which uniform and non-isotropic distributions have the same free energy) in all these 
numerical simulations. In particular, pi is close to p* (see Figure H]), so in most of 
the cases where both uniform and non-isotropic distribution are stable, the uniform 
distribution is not the global minimizer of the free energy, but in practical, meta- 
stability is very strong, and the solution still converges to the uniform distribution. 

5 Conclusion 

In this note, we have given the summary of strong results on the stability and 
instability of the equilibrium states of the modified Smoluchowski equation ([T]). 
This allows to have a precise description of the dynamics of the solution when time 
goes to infinity: it converges exponentially fast to a fixed equilibrium, with explicit 
formulas for the rates of convergence. We have also exhibited a specific example in 
which we observe a first order phase transition with a hysteresis loop (in contrast 
with the second order phase transition of the original Smoluchowski equation with 
dipolar potential [H]). The details of the proofs will be found in a longer paper [6] 
as well as numerical comparisons between the particle and kinetic models to confirm 
that the hysteresis is really intrinsic to the system and not simply an artifact of the 
kinetic modeling. 
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Figure 6: Hysteresis loop for the concentration parameter c in a numerical simulation 
of the kinetic equation ( IT6|) . with time varying p, in dimension 2. The red curve is 
the theoretical curve, the blue one corresponds to the simulation. 
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